Towards a Fully Ab-initio Description of the Diluted Magnetic Semiconductor 
Ga\- x Mn x As. Ferromagnetism, Electronic Structure, and Optical Response. 
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There are two competing views of itinerant ferromagnetism, the first viewing ferromagnetism as 
resulting from the indirect coupling between local moments via the itinerant carrier dynamics, 
the so-called RKKY mechanism, while in the alternative picture, ferromagnetism results from the 
spin polarization of itinerant carriers by the strong atomic Hund interaction - the so-called double 
exchange (DE) scenario. Which view describes the ferromagnetism in diluted magnetic semicon- 
ductors, materials with promise for spintronic applications, is still unclear. Here, we describe the 
detailed physical response of the prototype material Gai- x Mn x As using a combination of first- 
principles bandstructure with methods based on dynamical mean field theory to incorporate strong, 
dynamical correlations and intrinsic as well as extrinsic disorder in one single theoretical picture. 
We show how ferromagnetism is driven by DE, in agreement with very recent observations, along 
with a good quantitative description of the details of the electronic structure, as probed by scanning 
tunnelling microscopy (STM) and optical conductivity. Our results show how ferromagnetism can 
be driven by DE even in diluted magnetic semiconductors with small carrier concentration. 
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Spintronics is a rapidly emerging technology where it is 
not the electron charge but the electron spin that carries 
information. It offers promising opportunities for devel- 
oping a new generation of devices based on a combina- 
tion of standard microelectronics with spin-dependent ef- 
fects arising from the interaction of the carrier spin with 
the magnetism of the material [1]. Another promising 
area envisions its marriage with optical photons, with 
the possibility of designing new spin based devices such 
as spin- field effect transistor (FET), spin- light emitting 
diodes (LED), optical switches operating in the THz 
range, modulators, encoders and decoders, and bits for 
quantum computation [2] . Discovery of optically induced 
ferromagnetism in some Mn-doped III-V semiconductors 
opens up the possibility of applications to photonic stor- 
age devices and photonically-driven micromechanical el- 
ements [1]. Success of this vision demands a deeper un- 
derstanding of fundamental spin interactions in the solid 
state, along with the realistic bandstructures, and roles 
of dimensionality as well as intrinsic (extrinsic) defects. 

The discovery of dilute magnetic semiconductors 
(DMS), which are III-V semiconductors randomly doped 
with small amounts of magnetic atoms, like Mn 2+ , has 
provided us with an attractive example of prototype ma- 
terials of great interest in this context [3]. These have 
ferromagnetic transition temperatures much higher than 
those of earlier known E'u-chalcogenides [3], of order 
100 K, and the magnetism can be controlled electron- 
ically The optimum value of x in Ga\- x Mn x As cor- 
responding to highest T c = 100 K is 0.043 - 0.05. In 
spite of much activity, the nature of ferromagnetism, as 
well as the detailed electronic structure, is not properly 
understood. Ultimately, the mechanism of FM should 
go hand-in-hand with the details of spin interactions in a 
system of dilute Mn 2+ ions doped randomly in the GaAs 



host, and to the details of modification of the electronic 
structure due to Mn doping. 

Application of the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) mechanism provides some understanding of 
magnetic and transport properties of DMS [4], but is 
somewhat questionable in this case, since it is only valid 
when the exchange interaction, J « Ep, the Fermi en- 
ergy of the carriers. This is obviously not the case in 
DMS, where the carrier concentration is low [5] . An alter- 
native double-exchange (DE) model based approach [5] 
has recently been invoked within a model bandstructure 
to study the evolution of T c with x. However, the specific 
conditions under which DE ideas can be applied to DMS 
remain to be elucidated. 

Here, we discuss an implementation of DMFT which 
permits the incorporation of the actual electronic struc- 
ture of real materials into the DMFT equations. Using 
this method, we show how a detailed understanding of 
half-metallic ferromagnetism, along with a quantitative 
description of the correlated electronic structure of Mn- 
doped GaAs is achieved. 

GaAs is a well known band semiconductor. Its elec- 
tronic structure, however, is sensitive to As antisites 
forming in the bulk during growth. The importance of 
these defects in the GaAs host, constituting a source of 
intrinsic disorder, will become clear below. Each Mn 
ion in GaAs serves a dual purpose, acting as an ac- 
ceptor as well as a magnetic impurity. In reality, the 
situation is a bit more complex when the realistic lo- 
cal characteristics of the Mn ion, such as the charge 
transfer energy (Emu — Eas), an d the p — d hybridiza- 
tion (As — Mn), and Coulomb interactions, are taken 
into account. The importance of including such quan- 
tum chemical aspects is shown very clearly in spectro- 
scopic measurements [6] ; these yield the on-site Hubbard 
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U = 3.5 eV, the charge transfer energy A = —1.5 eV, 
and the hybridization t pc t = 1.1 eV. Additionally, the 
Hund's rule coupling Jh — 0.55 eV, acting like a mag- 
netic impurity potential [5]. Finally, the random distri- 
bution of the Mn ions gives rise to an additional non- 
magnetic disorder component. All previous theoretical 
works done for Ga\- x Mn x As have focused either on pure 
bandstructure aspects [10], or on model approaches [5], 
and we are not aware of calculations which incorporate 
the strong correlation at Mn sites, along with local quan- 
tum chemical information (t p( i, A) and intrinsic As anti- 
site disorder as well as Mn-doping related magnetic and 
non-magnetic disorder into the complex bandstructure of 
this material in a consistent way. Very recent theoreti- 
cal advances [11] have opened an attractive possibility to 
study the correlated bandstructure of materials involving 
transition metal oxides and rare-earth based compounds 
via LDA+DMFT. 

In what follows, we show how a consistent theoretical 
(via LDA+DMFT) implementation including real band- 
structure of GaAs, quantum chemical and strong correla- 
tion aspects mentioned above, and As antisite disorder is 
indeed necessary for a detailed quantitative understand- 
ing of the physics of Ga\- x Mn x As. 

A beautiful recent time-resolved magneto-optic mea- 
surement clearly shows the half-metallic character of the 
ferromagnetic metallic state in Ga\- x Mn x As [7]. But 
this is not all, as the correlated nature of this state is 
further revealed by various responses. First, photoemis- 
sion measurements [8] reveal very small spectral weight 
near the Fermi surface, and a T-dependent build-up of 
spectral weight at higher energies. Further, the change of 
the chemical potential Ep with Mn doping is intriguing. 
At small x, it appears to be almost pinned to its x — 
value, but starts moving towards the valence band (VB) 
around the x value where the semiconductor-metal tran- 
sition accompanied by FM occurs. At x > 0.05, it moves 
up again, concomitant with occurence of insulating be- 
havior with reduction of T c , clearly showing the intimate 
connection between changes in electronic structure and 
variation in T c as a function of x. More anomalous fea- 
tures are observed in optical measurements [9]: there is 
no quasicoherent Drude response; instead, the spectrum 
exhibits a peak centered around 1200 cm' 1 , with a curi- 
ous bump at lower energies. From a "modified Drude fit" , 
an effective mass m* = 0.72m is deduced, in contrast to 
nib = 0.24m from LDA calculations. The depression of 
the low-energy incoherent spectral weight with increasing 
T (decreasing magnetization) is also very clearly seen; 
what is interesting is that this spectral weight transfer 
(SWT) takes place over a wide energy scale of about 4 eV 
(much larger than T c ). These observations are inexpli- 
cable within bandstructure ideas, and constitute direct 
evidence for the importance of dynamical electronic cor- 
relations in Ga\- x Mn x As. 
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FIG. 1. Total density of states (DOS) of pure (without an- 
tisites) GaAs computed within the Kohn-Sham local density 
approximation [10]. The Fermi level is located at u = 0.0. 

While it provides a good description of ground state 
properties, the local density approximation (LDA) can- 
not describe the excited states of correlated systems, 
since the dynamical correlations in a quantum many- 
body system are neglected in the LDA. Thus, given the 
above arguments, a realistic description of Ga\- x Mn x As 
requires combining the LDA, which provides an excel- 
lent ab-initio one-electron bandstructure of weakly cor- 
related metals and band insulators [10,11], with well con- 
trolled, state-of-the-art many body calculations capable 
of accessing the dynamical effects of strong electron cor- 
relations in a proper way. Recently, the combination of 
LDA with dynamical mean field theory (DMFT) has been 
shown to provide good quantitative description for vari- 
ous correlated systems [12]. Using LDA+DMFT, we will 
show how good quantitative agreement with DE ferro- 
magnetism, STM data, and optical conductivity is ob- 
tained. 

Local atomic-like features of Mn 2+ in GaAs are ob- 
tained from the set of many-body levels describing pro- 
cesses of electron addition and/or removal in the atomic 
(i-shell. The constants U, Jh can be computed from first 
principles. Itinerant aspects of the actual solid enter via 
a dynamical bath function describing p — d hybridization; 
it quantifies the degree of itinerance of the d electrons, 
and needs to be determined in a selfconsistent way from 
the DMFT equations. With this, the problem is reduced 
to solving an (asymmetric, in general) multi-orbital An- 
derson impurity problem selfconsistently embedded in a 
dynamical bath, giving us the impurity Green's function 
Gimp{u) and the self-energy £(o>)- 

The bands of GaAs are obtained from the eigenvalues 
of the matrix one-electron Hamiltonian H(k), obtained 
from the LDA Kohn-Sham hamiltonian [11]. The actual 
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LDA density of states (DOS) for GaAs (without anti- 
sites) is shown in Fig. 1. 

To begin with, we model the effect of the random As 
antisite potential in pure GaAs using the usual coher- 
ent potential approximation (CPA) [13]. Next, given the 
small concentration of Mn sites, we first solve the asym- 
metric Anderson impurity problem in the dynamical bath 
provided by the GaAs bandstructure for a single Mn im- 
purity with the parameters given earlier, using the gener- 
alization of the iterated perturbation theory (IPT) for ar- 
bitrary filling [12]. The exact low- frequency behavior is 
obtained from the Friedel-Luttinger sum rule [12], while 
the correct high-frequency behavior is obtained from a 
selfconsistent computation of moments. These equations 
have appeared previously in various contexts [12], and we 
do not reproduce them here. For the one-band Hubbard 
model, very good agreement with quantum Monte Carlo 
data has been reported at high T . For a concentration 
x of Mn impurities, one has additionally to perform a 
configurational average over random Mn positions (non- 
magnetic) as well as an average over all spin configura- 
tions of the Mn ions (magnetic) . These effects of strong 
magnetic and positional disorder scattering are treated 
using the extended dynamical coherent potential approx- 
imation (CPA) [5], which is combined with IPT for the 
Mn impurities in a consistent way [14]. 

Marriage of LDA with DMFT gives us a quantitative 
description of both ground- and excited state properties. 
Finite temperature, and local moment effects above T c 
are readily handled using temperature Green's functions 
within DMFT. 

We now describe our results. In Fig. 2, we show the 
total density of states (DOS) for Ga\- x Mn x As for x = 
0, 0.022, 0.043 for our chosen parameter set. Given [15] 
that each Mn dopant creates an additional As anti-site, 
the concentration of antisites, n as — 0.015+a;. For x = 0, 
the effect of the random As anti-site disorder is modeled 
using CPA, resulting in an impurity (As) band of anti- 
sites split off from the VB in the semiconducting gap 
of GaAs. We draw attention to the good agreement of 
the calculated (x = 0) DOS with STM data [15]. For 
small x = 0.022, the total DOS still indicates semicon- 
ducting behavior, in coincedence with conclusions from 
Ref. [16]. Notice, however, the appearance of a second, 
broader peak in the gap. Again, this is in accordance 
with results of STM measurements [15]. Without ran- 
dom magnetic and non-magnetic impurity scattering, a 
quasicoherent peak is obtained (not shown). Inclusion of 
strong magnetic (Jh) and positional disorder scattering 
by CPA washes out this low-energy coherence, result- 
ing in the red curve in Fig. 2. For x — 0.043, metallic 
behavior, together with ferromagnetism (see below), is 
clearly revealed in the results (Fig. 2, green). From the 
real part of the one particle self-energy (not shown), we 
find 9S'(a;)/9w| w=/J = —2.03, leading to a mass enhance- 
ment of m*/mb = 3.03 (here, is the carrier mass 



estimated from LDA), in excellent agreement with ob- 
servations. The half-metallic character of the system for 
< 0.2 is clearly seen from our computed results (in- 
set to Fig. 2), in full accord [7] with recent experimental 
work. 
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FIG. 2. Renormalized total DOS of Gai- x Mn x As as a 
function of x. For x = 0, the additional impurity feature 
at uj — —0.1 is related to the antisite (n a3 = 0.015) contribu- 
tion (blue). For x — 0.022, a second Mn-doping related peak 
is resolved, but the system remains insulating (red). With 
x — 0.043 (green), half-metallic behavior is clearly seen in the 
inset (bold and dashed lines correspond to majority- and mi- 
nority-spin DOS), corresponding to DE ferromagnetism. No- 
tice the good agreement with STM data [15]. 

The ferromagnetic character of the metallic state is 
clearly shown by using the DMFT propagators to com- 
pute the static part of the spin susceptibility, Xs(q) = 
A -1 ^ k Gf(k + q)G|(k), and the effective intersite ex- 
change, J e //(q) = J 2 Xs(q)- Given the small concen- 
tration of the correlated Mn ions, we expect that ver- 
tex corrections appearing normally in the computation 
of Xs (q) will not modify our estimate of T c much (these 
corrections scale with the Mn concentration). The tran- 
sition temperature T c is estimated from the DMFT equa- 
tions [5]. We find, in full agreement with [16], that T c = 
for x = 0, 0.022 and T c = 120 K for x = 0.043. 

In light of our results, we see clearly that ferromag- 
netism in Ga\- x Mn x As results from the interplay be- 
tween two scales: the effective impurity bandwidth of 
the correlated solution (W e //) and the local Hund in- 
teraction (Jh) with Jh > We//, putting the system in 
the DE class. As described above, this is out of scope 
of pure LDA- or LDA+U based approaches, which can- 
not access dynamical effects of strong electronic correla- 
tions and static (dynamic) disorder. A proper treatment 
of these dynamical processes plays a crucial role (via 
DMFT) in generating an impurity band with reduced 
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bandwidth (W < Jh), stabilizing DE ferromagnetism in 
this DMS system. 
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FIG. 3. Computed optical conductivity a(ui,T) versus u> 
for various T both below and above T c — 120 K. Good quan- 
titative agreement with experiment [9] is clearly visible. 

Finally, we compute the optical conductivity, c(u>), 
directly from the fully renormalized DOS using the re- 
sult known rigorously [17] in DMFT. In Fig. 3, we show 
our result, a(u>,T), as a function of u> for x — 0.043. 
Concentrating on the T dependence, we see clearly that 
our results are very similar to the one from Ref. [9]. 
In particular, the broad peak around 800 cm -1 , the 
smaller (anti-site related) bump around 100 cmT 1 (some- 
what different from experiment), and the incoherent low- 
energy response are all in complete accordance with ob- 
servations. More satisfyingly, the T dependence of the 
spectral weight transfer is also correctly reproduced: at 
low-T (below T c ), no crossing point in <t(lu) is seen till 
3000 cm' 1 , while above T c , the curves seem to cross 
around 800 cm -1 , in full semiquantitative agreement 
with experiment. However, we were not able to resolve 
a clear crossing point in the spectra above T c . Finally, 
the distribution and T dependence of the optical spectral 
weight (decrease with increasing T below T c , and increas- 
ing with T above T c ) is in good agreement with the data 
as well. 

Given the detailed quantitative agreement of our re- 
sults with those gleaned from various spectroscopic, opti- 
cal, and magnetic measurements, we believe that we have 
provided the first totally ab-initio (LDA+correlations) 
description of Ga\- x Mn x As. Our study highlights the 
importance of including intrinsic and extrinsic disorder, 
strong dynamical correlations and magnetic scattering, 
in concert with the real LDA bandstructure in one single 
picture. We have clearly shown how a consistent treat- 
ment of dynamical electronic correlations and disorder 



generates a narrow, quasicoherent impurity band in the 
semiconductor band gap, opening up the possibility for 
DE to drive the system into a ferromagnetic half-metallic 
state. Other diluted magnetic semiconductors of great 
current interest can be studied readily within this frame- 
work. 
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